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Abstract In these lectures I discuss Feynman graphs and the associated Feynman 
integrals. Of particular interest are the classes functions, which appear in the eval- 
uation of Feynman integrals. The most prominent class of functions is given by 
multiple polylogarithms. The algebraic properties of multiple poly logarithms are 
reviewed in the second part of these lectures. The final part of these lectures is 
devoted to Feynman integrals, which cannot be expressed in terms of multiple poly- 
logarithms. Methods from algebraic geometry provide tools to tackle these integrals. 



1 Feynman graph polynomials 

The first part of these lectures is centred around two graph polynomials. We will 
give four different definitions of these two polynomials, each definition will shed a 
different light on the nature of these polynomials. The presentation in this section 
follows 0. 



1.1 Graphs 

Let us start with a few basic definitions: A graph consists of edges and vertices. We 
will mainly consider connected graphs. The valency of a vertex is the number of 
edges attached to it. Vertices of valency 0, 1 and 2 are special. A vertex of valency 
is necessarily disconnected from the rest of graph and therefore not relevant for 
connected graphs. A vertex of valency 1 has exactly one edge attached to it. This 
edge is called an external edge. All other edges are called internal edges. In the 
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physics community it is common practice not to draw a vertex of valency 1, but just 
the external edge. A vertex of valency 2 is called a mass insertion and is usually 
not considered. Therefore in physics it is usually implied that a genuine vertex has 
a valency of three or greater. 

An edge in a Feynman graph represents a propagating particle. The edges are 
drawn in a way as to represent the different types of particles. For example, one uses 
lines with an arrow for fermions, wavy lines for photons or curly lines for gluons. A 
simple line without decorations is used for scalar particles. To each (orientated) edge 
we associate a D-dimensional vector q and a number m, describing the momentum 
and the mass of the particle. D is the dimensions of space-time. 

Vertices of valency n > 3 represent interactions of n particles. At each vertex 
we have momentum conservation: The sum of all momenta flowing into the vertex 
equals the sum of all momenta flowing out of the vertex. 

To each Feynman graph we can associate a new graph, obtained by replacing 
each propagator of the original graph by a scalar propagator. This new graph is 
called the underlying topology. This new graph does no longer carry the informa- 
tion on the type of the particles propagating along the edges. We will later associate 
to each Feynman graph an integral, called the Feynman integral of this graph. It 
turns out, that the Feynman integral corresponding to an arbitrary Feynman graph 
can always be expressed as a linear combination of Feynman integrals correspond- 
ing to Feynman graphs with scalar propagators. Therefore it is sufficient to restrict 
ourselves to the underlying topology and to restrict our study to Feynman graphs 
with scalar propagators. 

Let us now consider a graph G with n edges and r vertices. Assume that the graph 
has k connected components. The loop number I is defined by 

l=n-r+k. (1) 

If the graph is connected we have l = n — r+\. The loop number / is also called the 
first Betti number of the graph or the cyclomatic number. In the physics context it 
has the following interpretation: If we fix all momenta of the external lines and if 
we impose momentum conservation at each vertex, then the loop number is equal 
to the number of independent momentum vectors not constrained by momentum 
conservation. 

A connected graph of loop number is called a tree. A graph of loop number 0, 
connected or not, is called a forest. If the forest hast k connected components, it is 
called a fc-forest. A tree is a 1 -forest. 



1.2 Spanning forests 

Given an arbitrary connected graph G, a spanning tree of G is a subgraph, which 
contains all the vertices of G and which is a tree. In a similar way, given an arbitrary 
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Fig. 1 An example of a Feynman graph and the associated set of spanning trees. 



connected graph G, a spanning &;-forest of G is a subgraph, which contains all the 
vertices of G and which is a A;-forest. 

We have already associated to each edge a momentum vector and a mass. In 
addition we associate now to each internal edge ej a real (or complex) variable xi. 
The variables xi are called Feynman parameters. For each graph we can define two 
polynomials "i/ and & in the variables Xj as follows: Let G be a connected graph 
and £7\ the set of its spanning trees. The first graph polynomial is defined by 

E U x j ( 2 ) 

This is best illustrated with an example. Fig. ([D shows a Feynman graph decorated 
with the Feynman parameters x\ to X5, as well as the associated set 3F\ of span- 
ning trees. For each spanning tree we take the Feynman parameters associated to 
the edges not belonging to this spanning tree. Summing over all spanning trees we 
obtain for this example 

% = X\X2 + + JC1JC3 + X2X4 + X2X5 + X\X$ + X4X5 + . (3) 

is also called the first Symanzik polynomial of the graph G. In mathematics, the 
Kirchhoff polynomial of a graph is better known. It is defined by 

X= £ Y\ Xj . (4) 

The difference between the two definitions is given by the fact that in the case of 
J(T we consider all edges belonging to the spanning tree T, while in the case of tyt 
we consider all edges not belonging to T. There is a simple relation between the 
Kirchhoff polynomial Jtf and the first Symanzik polynomial <fy : 

W(x u ...,x n ) =x\...x n X (—,...,— ) . (5) 

We now turn to the definition of JF. Let G be a connected graph and 3*2 me set of 
its spanning 2-forests. An element of 5^ is denoted as (T\ , T2). Let us further denote 
by Pjj the set of external momenta of G attached to 7}. We first define a polynomial 

=^oby 
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Fig. 2 An example of a Feynman graph and the associated set of spanning 2-forests. 



(Ti,T 2 )eS 2 



n * 



Pj-Pk 



(6) 



Here, /? ; - • is the Minkowski scalar product of two momenta vectors, fi is an arbi- 
trary scale introduced to make the expression dimensionless. J£" is defined by 



2 



(7) 



OTj denotes the mass of the ;-th internal line. If all internal masses are zero, we have 
& = ,j£q. c£ is called the second Symanzik polynomial. Again, let us illustrate the 
definition of with an example. Fig. (|2]i shows a Feynman graph and the associ- 
ated set S?i of 2-forests. For simplicity we assume that all internal masses are zero, 
therefore we have & = J^o- In fig. (f2]i we have labelled the internal edges with the 
Feynman parameters x\ to X4, The external edges have been labelled with the exter- 
nal momenta p\ to p^. We orientate these edges such that p\ to P4 are all flowing 
outwards. With this choice momentum conservation reads 



Pi +P2 + P3+P4 = 0. 

The Mandelstam variables s and t are defined by 



*= (P1+P2) 2 , t = (p 2 +p 3 ) 2 . 



(8) 



(9) 



From the definition in eq. © we find for the polynomial & for this example 



(-*) (-0 (-pi) , . . (-pi) , . „ (-pi) , „ . (-pI) 



jLi 2 



X1X4- 



-X\X2- 



-X2X3- 



-X3X4- 



A few remarks on the two Symanzik polynomials are in order: Both polynomials 
are homogeneous in the Feynman parameters, tyt is of degree I, is of degree 
/ + 1 . The polynomial is linear in each Feynman parameter. If all internal masses 
are zero, then also & is linear in each Feynman parameter. In expanded form each 
monomial of % has coefficient +1. 
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1.3 Feynman integrals 



Feynman graphs have been invented as a pictorial notation for mathematical ex- 
pressions arising in the context of perturbative quantum field theory. Each part in 
a Feynman graph corresponds to a specific expression and the full Feynman graph 
corresponds to the product of these expressions. For scalar theories the correspon- 
dence is as follows: An internal edge corresponds to a propagator 

~2~ — 2 ' ( 10 > 

q L — m L 

an external edge to the factor 1. A vertex corresponds in scalar theories also to 
the factor 1. In addition, there is for each internal momentum not constrained by 
momentum conservation an integration 

d D k 

Let us now consider a Feynman graph G with m external edges, n internal edges 
and / loops. To each internal edge we associate apart from its momentum and its 
mass a positive integer number v, giving the power to which the propagator occurs. 
(We can think of v as the relict of neglecting vertices of valency 2. A number v > 1 
corresponds to v — 1 mass insertions on this edge). The momenta flowing through 
the internal lines can be expressed through the independent loop momenta k\, k[ 
and the external momenta p\, p m as 

Z m 

Y.l ): ' k > ■ E°W»7j P<h a <i e {-1,0,1}. (12) 

j=l 7=1 

We define the Feynman integral by 

n r ( V ^ / jDi n 

*>- rir^j <"> /n-?n^ <13) 

with v = Vi +... + V„. The prefactorin front of the integral is the convention used in 
this article. This choice is motivated by the fact that after Feynman parametrisation 
we obtain a simple formula. Feynman parametrisation makes use of the identity 

j=' p j nr(v,)i 




where co is a differential (n — l)-form given by 
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C0=52(-l)-' 1 Xj dxi A ... A dxj A ... A dx„ . (15) 

The hat indicates that the corresponding term is omitted. The integration is over 

A = {[xi :x 2 : ...:x„] e P"" 1 |jt; > 0, 1 </<«}. (16) 
We use eq. ( TBI ) with Pj = —q 1 - + rrpj. We can write 

n ill 

E Xj(-qj + m)) EE +E2^ r .g r - /, (17) 

/=1 r=l.v=l r=l 

where M is a / x / matrix with scalar entries and Q is a /-vector with D-vectors as 
entries. After Feynman parametrisation the integrals over the loop momenta k\, 
k/ can be done and we obtain 

a/v-(l+l)D/2 

(18) 

The functions and & are given by 

f/=det(M), ^ = det(M)(-y + gM" 1 e)/ J u 2 . (19) 

It can be shown that eq. < fT~9b agrees with the definition of ^ and & given in sec- 
tion \L2\ in terms of spanning trees and spanning forests. Thus, eq. dT9b provides a 
second definition of the two graph polynomials. Eq. ([T81 defines the Feynman inte- 
gral of a graph G in terms of the two graph polynomials & and A few remarks 
are in order: The integral over the Feynman parameters is a (n — 1) -dimensional in- 
tegral in projective space P"~ 1 , where n is the number of internal edges of the graph. 
Singularities may arise if the zero sets of % and & intersect the region of integra- 
tion. The dimension D of space-time enters only in the exponents of the integrand 
and the exponents act as a regularisation. 




1.4 The Laplacian of a graph 

For a graph G with n edges and r vertices define the Laplacian L l42l[38l as a r x r- 
matrix with 

_ ( Y^ x k if i = j an d edge is attached to V; and is not a self-loop, 
lJ \— Y,*k if ! 7^ j an d edge connects v, and vj. 

We speak of a self-loop (or tadpole) if an edge starts and ends at the same vertex. 
In the sequel we will need minors of the matrix L and it is convenient to introduce 
the following notation: For a r x r matrix A we denote by A[j'i, ...,ik;ji , A] tne 
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(r — k) x [r — k) matrix, which is obtained from A by deleting the rows i\, and 
the columns j\, j^. For A[/i,...,4;m, •■•,'*] we w iU simply write A[ii, ...,ijt]. The 
matrix-tree theorem relates the Laplacian of a graph to its Kirchhoff polynomial: 

jT = detL[z]. (20) 

A generalisation by the all-minor matrix tree theorem lfP3] H4l |3~T1 leads to the fol- 
lowing expressions for the graph polynomials & and JPq: Starting from a graph G 
with « internal edges, r internal vertices (vi , v r ) and m external legs, we first at- 
tach m additional vertices (v r+ i, ...,v r + m ) to the ends of the external legs and then 
associate the parameters z\, z m with the external edges. This defines a new graph 
G. We now consider the Laplacian L of G and the polynomial 

W(xi,...,x n ,zi,...,z m ) = det L[r+l,...,r + m]. (21) 

We then expand W in polynomials homogeneous in the variables zf 

w = W {0) + W W +w {2) + ... + w { - m \ 

W { k ) = ^ ^ Ji) (x l ,...,x n )zj 1 ...z Jk . (22) 

l<jl<-<j k <m 

We then have 

= W {0 \ 

&=x 1 ..jc n W$ f^-j-.^-J for any;, 

This provides a third definition of the Feynman graph polynomials °i/ and This 
formulation is particularly well suited for computer algebra. 



1.5 Deletion and contraction properties 

Let us now consider a recursive definition of the two graph polynomials based on 
deletion and contraction properties. We first define a regular edge to be an edge, 
which is neither a self-loop nor a bridge. In graph theory an edge is called a bridge, 
if the deletion of the edge increases the number of connected components. Examples 
for graphs containing either a self-loop or a bridge are shown in fig. (O. For a graph 
G and a regular edge e we define 

G/e to be the graph obtained from G by contracting the regular edge e, 
G — e to be the graph obtained from G by deleting the regular edge e. (24) 
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Fig. 3 Examples of graphs 
containing a self-loop (left) or 
a bridge (right). 





self-loop 



bridge 



Fig. 4 A graph G, together 
with the graph G — e, where 
the edge e has been deleted 
and the graph G/e, where the 
edge e has been contracted. 




The operations of deletion and contraction are illustrated in fig. ©. For any regular 
edge e k we have 

<V(G) = W(G/e k )+x k W(G-e k ), 

& (G) = & (G/e k )+x k & (G-e k ). (25) 

The recursion terminates when all edges are either bridges or self-loops. These 
graphs are called terminal forms. If a terminal form has r vertices and / (self-) loops, 
then there are (r — 1) "tree-like" propagators, where the momenta flowing through 
these propagators are linear combinations of the external momenta pi alone and in- 
dependent of the independent loop momenta kj. The momenta of the remaining / 
propagators are on the other hand independent of the external momenta and can be 
taken as the independent loop momenta kj, j = I,..., I. Let us agree that we label the 
{r — 1) "tree-like" edges from 1 to r— 1, and the remaining / edges by r, n with 
n = r + 1 — 1. We further denote the momentum squared flowing through edge j by 
q 1 -. For a terminal form we have 

r-i ( -q 2 \ 

<?/ =x r ...x n , J^o =x r ...x„Y,Xj — y ■ ( 26 ) 



7=1 



Eq. d25l l together with eq. (|26*| | provides a fourth definition of the graph polynomials 
9/ and 

Let e a and e b be two regular edges, which share a common vertex. We have the 
following factorisation theorems: 

(G/e a - e b ) W (G/e b - e a ) -W{G-e a - e b ) W (G/e a /e b ) = 



W (G/e a - e b ) & (G/e b - e a ) -<%{G- e a - e b ) J*b (G/e a /e b ) 

(G/e a - e b ) "7/ (G/e b - e a ) - J*b (G - e a - e b ) W (G/e a /e b ) = 
(Ai_\ ( A 2 



2 — — ] . (27) 
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A\ and A2 are polynomials in the Feynman parameters and can be expressed as sums 
over 2-forests and sums over 3-forests, respectively []8]. If for all external momenta 
one has 

(j>h ■ Ph.) ■ (Ph ■ Pk) = (Ph • Pis) • (Ph ■ Pk) > ( 28 > 

then 



(G/e a - e b ) &o (G/e b -e a )- & (G - e a - e b ) ^ (G/e a /e b ) = 




The factorisation theorems follow from Dodgson's identity [15. 48], which states 
that for any n x n matrix A one has 

det(A)det(A[/,;]) = det (A[i]) det (A[j}) - det (A [/; y]) det {A\J;i]) . (29) 

We recall that 

A [i] is obtained from A by deleting the z-th row and column, 
A [i; j) is obtained from A by deleting the ;-th row and the y'-th column, 
A[i,j] is obtained from A by deleting the rows and columns ; and j. 

The first formula of eq. (l27l is at the heart of the reduction algorithm of lfTT1[T2l . 



2 Multiple polylogarithms 

Let us come back to the Feynman integrals defined in eq. (fT8l . A Feynman integral 
has an expansion as a Laurent series in the parameter e = (4 — D)/2 of dimensional 
regular! sation: 

I G = £ c j£ J. (30) 

;=-2/ 

The Laurent series of an /-loop integral can have poles in e up to the order (2/). 
The poles in e correspond to ultraviolet or infrared divergences. The coefficients cj 
are functions of the scalar products pj ■ p^, the masses m; and (in a trivial way) of 
the arbitrary scale fi. An interesting question is, which functions do occur in the 
coefficients Cj. 
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The question, which functions occur in the coefficients c\ has a satisfactory answer 
for one-loop integrals. If we restrict our attention to the coefficients cj with j < 
(i.e. to c-2, c_i and cq), then these coefficients can be expressed as a sum of 
algebraic functions of the scalar products and the masses times two transcendental 
functions, whose arguments are again algebraic functions of the scalar products and 
the masses. 

The two transcendental functions are the logarithm and the dilogarithm: 




Ii 2 (jc) = £ (31) 



2.2 The sum representation of multiple polylogarithms 

Beyond one-loop an answer to the above question is not yet known. We know how- 
ever that the following generalisations occur: From eq. d3H it is not too hard to 
imagine that the generalisation includes the classical polylogarithms defined by 

Li m (x) = £ £ (32) 

However, explicit calculations at two-loops and beyond show that a wider generali- 
sation towards functions of several variables is needed and one arrives at the multiple 
polylogarithms defined by Il22ll23l l9l 

°° x" [ x" k 

Li '"i m k (x u ...,x k )= £ * (33) 

«l>n 2 >— >«jt>0 1 * 

Methods for the numerical evaluation of multiple polylogarithms can be found in 
l44ll . The values of the multiple polylogarithms at x\ = ...Xk = 1 are called multiple 
C-values J9]|6): 

Cmi mj = Li mi: ,„ 2 ,...,m jt (l, 1, 1) = ^ ~ mf ' ••• ' ~ m£ ■ (34) 

«l>«2>- ->"*>0 "l n k 

Important specialisations of multiple polylogarithms are the harmonic polyloga- 
rithms 113611211 

Jt-1 
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Further specialisations leads to Nielsen's generalised polylogarithms 11341 

S n ,p(x) = Li„+i,i i(x,l,...,l). (36) 

p-i 

2.3 The integral representation of multiple polylogarithms 

In eq. d33l l we have defined multiple polylogarithms through the sum representation. 
In addition, multiple polylogarithms have an integral representation. To discuss the 
integral representation it is convenient to introduce for Zky^O the following functions 



tl 'k-l 
dt\ f dt2 f dtk 



h ~Z\ J t 2 ~Z2 J t k ~Zk 

oo o 



(37) 



In this definition one variable is redundant due to the following scaling relation: 

G{zi,...,Zk;y) = G(xzi,...,xzk;xy) (38) 
If one further defines g(z',y) = \/{y — z), then one has 

—G(zi,...,zkiy) = g(zi;y)G(z 2 ,-,Zk;y) (39) 

and 



G(zi,z2,-,z k ;y) = J dt g (zi;t)G(z2,...,z k ;t). (40) 
o 

One can slightly enlarge the set and define G(0, ...,0;;y) with k zeros for z\ to Zk to 
be 

G(0,...,0;y) = ^(\ny) k . (41) 

This permits us to allow trailing zeros in the sequence (zi,...,Zjt) by defining the 
function G with trailing zeros via eq. (l40t and eq. (HTt . To relate the multiple poly- 
logarithms to the functions G it is convenient to introduce the following short-hand 
notation: 

G mi m k {zu—,Zk\y) = G(Q^Q,zi,...,Zk-i,O^Q,Zk-,y) (42) 

m ] — 1 — 1 

Here, all Zj for j = 1, are assumed to be non-zero. One then finds 



12 Stefan Weinzierl 



Li m] m k {xi,-,x k ) = (-l)*G mi m4 ( — , — !— — - — ;1 ) . (43) 

* X\ x\x 2 X\...Xh 



The inverse formula reads 



G m , mt (zi,...,z it ;y) = (-l)*Li m] ( Z,2.,...,*_L ) . (44) 

Eq. d43l together with eq. (J42]> and eq. ( f37T > defines an integral representation for the 
multiple polylogarithms. As an example, we obtain from eq. (|43l and eq. d38l ) the 
integral representation of harmonic polylogarithms: 

H mi ,...,m k (x) = (-1) G mu ... m {l,...,l;x). (45) 

The function G m| mjt (l,..., l;x) is an iterated integral in which only the two one- 
forms 

®o = —, ®i = — (46) 

corresponding to z = and z = 1 appear. If one restricts the possible values of z 
to zero and the n-th roots of unity one arrives at the class of cyclomatic harmonic 
polylogarithms (TJ. 



2.4 Shuffle and quasi-shujfle algebras 

Multiple polylogarithms have a rich algebraic structure. The representations as it- 
erated integrals and nested sums induce a shuffle algebra and a quasi-shuffle alge- 
bra, respectively. Shuffle and quasi-shuffle algebras are Hopf algebras. Note that the 
shuffle algebra of multiple polylogarithms is distinct from the quasi-shuffle algebra 
of multiple polylogarithms. 

Consider a set of letters A. The set A is called the alphabet. A word is an ordered 
sequence of letters: 

w = hl 2 ...l k . (47) 

The word of length zero is denoted by e. Let K be a field and consider the vector 
space of words over K. A shuffle algebra on the vector space of words is defined 
by 

{hh—k) ■ (k+l—lr) = 2- ^<7(l) ; cr(2)---'<7(r), ( 4 8) 



where the sum runs over all permutations a, which preserve the relative order of 
1,2,..., k and of lc+ l,...,r. The name "shuffle algebra" is related to the analogy 
of shuffling cards: If a deck of cards is split into two parts and then shuffled, the 
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relative order within the two individual parts is conserved. A shuffle algebra is also 
known under the name "mould symmetral" lH6l . The empty word e is the unit in 
this algebra: 

e ■ w = w ■ e = w. (49) 
A recursive definition of the shuffle product is given by 

(hh-k) ■ (k+l-lr) = h [(k-k) ■ (k+l-lr)] + k+1 [{hh-h) ■ (k+2-lr)] ■ (50) 

It is a well known fact that the shuffle algebra is actually a (non-cocommutative) 
Hopf algebra 1371 . In a Hopf algebra we have in addition to the multiplication and 
the unit a counit, a comultiplication and a antipode. The unit in an algebra can 
be viewed as a map from K to A and multiplication in an algebra can be viewed 
as a map from the tensor product A (g>A to A (e.g. one takes two elements from 
A, multiplies them and gets one element out). The counit is a map from A to K, 
whereas comultiplication is a map from A to A <g> A. We will always assume that the 
comultiplication is coassociative. The general form of the coproduct is 

A{a)=^ ] ®af\ (51) 
i 

where a,- denotes an element of A appearing in the first slot of A ®A and af^ 
correspondingly denotes an element of A appearing in the second slot. Sweedler's 
notation fl39l consists in dropping the dummy index i and the summation symbol: 

A (a) = fl (1) ®a (2) (52) 

The sum is implicitly understood. This is similar to Einstein's summation conven- 
tion, except that the dummy summation index i is also dropped. The superscripts ^ 
and ( 2 ) indicate that a sum is involved. Using Sweedler's notation, the compatibility 
between the multiplication and comultiplication is expressed as 

A(a-b)= («(')• b^)®(aW-bW). (53) 

The antipode S is map from A to A, which fulfils 

fl« ■ S (fl< 2 )) = S (a«) • a® = e ■ e(a). (54) 

With this background at hand we can now state the coproduct, the counit and the 
antipode for the shuffle algebra: The counit e is given by: 

«(«) = !, e{hl 2 ...l„)=0. (55) 



The coproduct A is given by: 
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Fig. 5 A shuffle algebra 
follows from replacing the 
integral over the square by 
an integral over the lower 
triangle and an integral over 
the upper triangle. 
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+ 



k 

A{hl 2 ...l k ) = £ (l J+1 ...l k )®(h...lj). (56) 

7=0 

The antipode S is given by: 

S(hl 2 ...l k ) = {-l) k hh-x-hh. (57) 

The shuffle algebra is generated by the Lyndon words. If one introduces a lexico- 
graphic ordering on the letters of the alphabet A, a Lyndon word is defined by the 
property w < v for any sub-words u and v such that w = uv. 

An important example for a shuffle algebra are iterated integrals. Let [a,b] be a 
segment of the real line and f\,f%, ... functions on this interval. Let us define the 
following iterated integrals: 

b t\ tk-1 

I(fi,fi,...,f k ;a,b) = J dhh{h) J dt 2 f 2 {t 2 )... J dt k f k (t k ) (58) 

a a a 

For fixed a and b we have a shuffle algebra: 

I(fi,f2,-,fkia,b)-I(f k+ i,...,f r ;a,b)= £ I(f a{l) ,f a{2) ,...,f a{r y,a,b), (59) 

shuffles a 

where the sum runs over all permutations a, which preserve the relative order of 
1,2, ...,k and of k+ l,...,r. The proof is sketched in fig. [5] The two outermost inte- 
grations are recursively replaced by integrations over the upper and lower triangle. 
The definition of multiple polylogarithms in eq. (f37T > is of the form of iterated inte- 
grals as in eq. d58l . Therefore it follows that multiple polylogarithms obey a shuffle 
algebra. An example for the multiplication is given by 

G(z\;y)G(z 2 ;y) = G(z\,z 2 ;y) + G{z 2 ,z\;y) ■ (60) 

Let us now turn to quasi-shuffle algebras. Assume that for the set of letters we have 
an additional operation 

(.,.) :A®A ->A, 

h®l 2 ->(h,h), (61) 



which is commutative and associative. Then we can define a new product of words 
recursively through 
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(hh-k) * (4+i -lr) = h [(h-h) * ik+i-lr)} + h+\ [(hh-h) * (k+i-l,-)} 

+ (l u l k+l )[(l 2 ...l k )*(h+2-lr)}- (62) 

This product is a generalisation of the shuffle product and differs from the recur- 
sive definition of the shuffle product in eq. d50t through the extra term in the last 
line. This modified product is known under the names quasi-shuffie product 11251 . 
mixable shuffle product [f24|. sniffle product or mould Symmetrel [[TBI . Quasi- 
shuffie algebras are Hopf algebras. Comultiplication and counit are defined as for 
the shuffle algebras. The counit e is given by: 

«(«) = !, e{hh-ln)=0. (63) 
The coproduct A is given by: 

k 

A{hl 2 ...l k ) = Y (/, ,.. ./a) •: /|. ../;;. (64) 

y'=o 

The antipode S is recursively defined through 

jfe-i 

S(hh-h) - -hl2...k-Z S ( l J+i- l k)*(h-lj), S(e) = e. (65) 

An example for a quasi-shuffie algebra are nested sums. Let n a and rif, be integers 
with n a < rib an d let /i> /2, ■■■ be functions defined on the integers. We consider the 
following nested sums: 

n b 'A — 1 — 1 

S(fi,f2,-,fk\n a ,n b ) = fd'i) Y, Mh)- Y /*(»*)• (66) 

i\=n a h=n a >k="a 

(The letter S denotes here a function, and not the antipode.) For fixed n a and iij, we 
have a quasi-shuffie algebra: 

S(fi,f2,-,fk;n a ,n b )*S(f k+1 ,...,f r ;n a ,n b ) = 

Y Mh) S(f 2 ,...,fk;n a ,h - \)*S(fk+i,—Jr\n a ,h - 1) 

h="a 
"b 

+ Y fkUi) s (fuf2,-,fk;n a ,ji - l)*S(f k + 2 ,...,f r ;n a ,ji - 1) 

j\="a 

"b 

+ Y /i(0/*(0 s (f2,—,fk;n a ,i-l)*S(fk +2 ,...,f r ;n a ,i-l) (67) 

i=n a 

Note that the product of two letters corresponds to the point-wise product of the two 
functions: 



(fiJj) (n) = fi(n)fj(n). 



(68) 
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• • • • 



Fig. 6 A quasi-shuffle algebra follows from replacing the sum over the square by a sum over the 
lower triangle, a sum over the upper triangle and a sum over the diagonal. 



The proof that nested sums obey the quasi-shuffle algebra is sketched in fig. ©. The 
outermost sums of the nested sums on the l.h.s of (167} are split into the three regions 
indicated in fig. (|6). The definition of multiple polylogarithms in eq. 03 is of the 
form of nested sums as in eq. ( 166} . Therefore it follows that multiple polylogarithms 
obey also a quasi-shuffle algebra. An example for the quasi-shuffle multiplication is 
given by 

Li m| (xj )Li,„ 2 (x 2 ) = Li mi ,„ 2 [x\ ,X2) + Lim2,mi (x2,x\ ) + Li m| _|_ m2 (jci*2). (69) 



2.5 Mellin-Barnes transformation 



In sect. ll.3l we saw that the Feynman parameter integrals depend on two graph poly- 
nomials °M and j£\ which are homogeneous functions of the Feynman parameters. 
In this section we will continue the discussion how these integrals can be performed 
and exchanged against a (multiple) sum over residues. The case, where the two 
polynomials are absent is particular simple: 

f\r(vj) 

r- (70) 




r(vi + ... + V n ) 

With the help of the Mellin-Barnes transformation we now reduce the general case 
to eq. (l70b . The Mellin-Barnes transformation reads 



1 1 



(Ai+A 2 + ...+A„) C ' = -F7ir ,„ i / doi... I dOn-i (71) 



xr(-cT 1 )...r(-cr„_ 1 )r( t T 1 + ... + tr„_ 1 + c)A^...A^r l l A- ai -- a '- l - c . 

Each contour is such that the poles of F(— a) are to the right and the poles of F(<7 + 
c) are to the left. This transformation can be used to convert the sum of monomials of 
the polynomials % and & into a product, such that all Feynman parameter integrals 
are of the form of eq. (l70t . As this transformation converts sums into products it is 
the "inverse" of Feynman parametrisation. With the help of eq. d70b we may perform 
the integration over the Feynman parameters. A single contour integral is then of the 
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form 




da 



r(a + a 1 )...r(a + a m )r(-a + bi)...r(-a + b„) 
r{a + c 2 )...r(a + c p ) r(-a + di)...r(-a + d q ) 



X 



■ (72) 



The contour is such that the poles ofF(<7 + fli), r(a + a m ) are to the right of 
the contour, whereas the poles of r(— a + b\), F(— a + b„) are to the left of the 
contour. We define 



Then the integral eq. (l72l converges absolutely for a > ifTTl and defines an analytic 
function in 



The integral eq. d72l i is most conveniently evaluated with the help of the residuum 
theorem by closing the contour to the left or to the right. Therefore we need to know 
under which conditions the semi-circle at infinity used to close the contour gives a 
vanishing contribution. This is obviously the case for |x| < 1 if we close the contour 
to the left, and for \x\ > 1, if we close the contour to the right. The case |jc| = 1 
deserves some special attention. One can show that in the case j3 = the semi-circle 
gives a vanishing contribution, provided X < — 1 . To sum up all residues which lie 
inside the contour it is useful to know the residues of the Gamma function: 

(-1)" (-1)" 
res (r(a + a),a = —a — n) = - — — , res (r(—a + a), a = a + n) = — - — — . 

nl «! 

In the general case, the multiple integrals in eq. (TTTb lead to multiple sums over 
residues. 



2.6 Z-sums 

The multiple sums over the residues can be expanded as a Laurent series in the 
dimensional regularisation parameter e. For particular integrals the coefficients of 
the Laurent series can be expressed in terms of multiple poly logarithms. To see this, 
we first introduce a special form of nested sums, called Z-sums I30, l45ll46ll29l : 



a = m + n — p — q, J3 =m—n — p + q, 




(73) 




(74) 



Z(n;m u ...,m k ;x u ...,x k ) 



f'l>!2>...>ijfc>0 



E 



n 




(75) 
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k is called the depth of the Z-sum and w — m\ + ... + m k is called the weight. If the 
sums go to infinity (« = °°) the Z-sums are multiple polylogarithms: 

Z(°°;m 1 ,...,m k ;x l ,...,x k ) = U mu ... >mk (xi,...,x k ). (76) 

For x\ = . . . = x k = 1 the definition reduces to the Euler-Zagier sums [[18] |47] |43] [7] 
®: 

Z(n;m\,...,m k ;\,...,\) = Z nn ,.,., mk (n). (77) 

For n = oa and x\ = ... = = 1 the sum is a multiple £-value: 

Z(oo;mi,...,m k ;l,...,l) = C>n m k - (78) 

The Z-sums are of the form as in eq. (lo*6l and form therefore a quasi- shuffle algebra. 
The usefulness of the Z-sums lies in the fact, that they interpolate between multiple 
polylogarithms and Euler-Zagier sums. Euler-Zagier sums appear in the expansion 
of the Gamma-function: 

r(n + e)=r(l + e)r(n) (79) 
x [l+eZ 1 (n-l) + e 2 Zi 1 («-l) + e 3 Z 1 n(n-l) + ... + e"- 1 Z 11 ...i(n-l)] . 

The quasi-shuffle product can be used to reduce any product 

Z(n;mi, ...;x\, ...) -Z(n\m\,...\x\,...) (80) 

of Z-sums with the same upper summation index n to a linear combination of single 
Z-sums. The Hopf algebra of Z-sums has additional structures if we allow expres- 
sions of the form 1301 

x n 

-f±-Z(n;m l ,...,m k ;x l ,...,x k ), (81) 
n o 

e.g. Z-sums multiplied by a letter. Then the following convolution product 

£ -Z(i-l;...)j^— 7 Z(n-i-V,...) (82) 
r-[ i m [n — i) m 

can again be expressed in terms of expressions of the form (1811 . In addition there is 
a conjugation, e.g. sums of the form 

-t l (■)(- 1 ) , ^•■•) (^) 

can also be reduced to terms of the form (18TV The name conjugation stems from 
the following fact: To any function f(n) of an integer variable n one can define a 
conjugated function Co /(«) as the following sum 
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Cof{n) = H n \{-\)'f(i). (84) 



= 1 

Then conjugation satisfies the following two properties: 

Col = 1, 

CoCof(n)=f(n). (85) 

Finally there is the combination of conjugation and convolution, e.g. sums of the 
form 

'7—1 / „ \ v i y n ~i 



-g^J(- 1 )'^---)^ S? ^-ft--.) (86) 

can also be reduced to terms of the form ( |8TV 

With the help of these algorithms it is possible to prove that the Laurent series in 
e of specific Feynman integrals contains only multiple polylogarithms (0. 



3 Beyond multiple polylogarithms 

Although multiple polylogarithms form an important class of functions, which ap- 
pear in the evaluation of Feynman integrals, it is known from explicit calculations 
that starting from two-loop integrals with massive particles one encounters func- 
tions beyond the class of multiple polylogarithms. The simplest example is given 
by the two-loop sunset diagram with non-zero masses, where elliptic integrals make 
their appearance ll28l . Differential equations provide a tool to get a handle on these 
integrals ||27l |26l l35l [T9l |2T1 |20l |2l [331 . 



3.1 The two-loop sunset integral with non-zero masses 

In this subsection we review how ideas of algebraic geometry can be used to ob- 
tain a differential equation for the Feynman integral ll32l . We first focus on the 
example of the two-loop sunrise integral. The two-loop sunrise integral is given in 
D-dimensional Minkowski space by 

, tonV (M 2 ) 3 '° fiPkuPt* 1 

^"'P > r(3-D)J in q in q {-k\ + m\) (-*f+m*) + ^ } 

with £3 = p — k\ — k2- Here we suppressed on the l.h.s. the dependence on the inter- 
nal masses m\, mi and and on the arbitrary scale ji. It is convenient to denote the 
momentum squared by t = p 2 . In terms of Feynman parameters the integral reads 
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Fig. 7 The intersection of 
the domain of integration A 
with the zero set X of the 
second Symanzik polynomial 
(left) and the elliptic curve 
y 2 — x 3 —x + I (right). 




r ^ 3 -i D 

S(D,t) = J a>-p5=D> ^ 

A 

where the two Feynman graph polynomials are given by 

= [— x\X2Xj,t + (x\m\ +X2^2 +xjm^j fy] /i~ 2 , % =x\x%- s rX2Xj,-\-X'iX\. 

It is simpler to consider this integral first in D = 2 dimensions and to obtain the result 
in D = 4 — 2e dimensions with the help of dimensional recurrence relations ll40l 
|4T1 . In two dimensions this integral is finite, depends only on the second Symanzik 
polynomial & and is given by 

S(2,t) =/J, (89) 

A 

From the point of view of algebraic geometry there are two objects of interest in 
eq. ( |89l : On the one hand the domain of integration A and on the other hand the 
algebraic variety X defined by the zero set of J? = 0. The two objects X and A 
intersect at the three points [1 : : 0], [0 : 1 : 0] and [0:0: 1] of the projective 
space P 2 . This is shown in fig. © on the left. We blow-up P 2 in these three points 
and we denote the blow-up by P. We further denote the strict transform of X by 
Y and the total transform of the set {X1X2X3 = 0} by B. With these notations we 
can now consider the mixed Hodge structure (or the motive) given by the relative 
cohomology group |0] 

H 2 (P\Y,B\BDY). (90) 

In the case of the two-loop sunrise integral considered here essential information on 
H 2 (P\Y,B\B n Y) is already given by H 1 (X). We recall that the algebraic variety X 
is defined by the second Symanzik polynomial: 

— xiX2X3f + (xim 2 +X2OT2 +X3m 2 ) (xi^2 +X2X3 +X3X1) =0. (91) 

This defines for generic values of the parameters t, m\, 1112 and an elliptic curve. 
The elliptic curve varies smoothly with the parameters t, mi, ni2 and 1113. By a bira- 
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tional change of coordinates this equation can brought into the Weierstrass normal 
form 

y 2 Z-x 3 -a 2 {t)xz 2 -a 3 (t)z 3 = 0. (92) 

The dependence of a 2 and a 3 on the masses is not written explicitly. In the chart 
z = 1 this reduces to 

y 2 -x 3 -a 2 (t)x-a 3 (t) = 0. (93) 

The curve varies with the parameter t . An example of an elliptic curve is shown in 
fig. (0 on the right. It is well-known that in the coordinates of eq. ( 1931 the cohomol- 
ogy group H 1 (X) is generated by 

dx d 
77 = — and 77 = — 77. (94) 
y at 

Since H l (X) is two-dimensional it follows that r) = ^rj must be a linear combina- 
tion of 77 and 77 . In other words we must have a relation of the form 

Po{t)Ti+pi(t)f}+p 2 (t)ri =0. (95) 
The coefficients po(t), p\{t) and p 2 (t) define the Picard-Fuchs operator 

LP) =p Q (t)^+p x (t)j t+P2 (t). (96) 
Applying the Picard-Fuchs operator to our integrand gives an exact form: 

i (2) (|)=^. (97) 

The integration over A yields 

L^S(2,t) = jdp = jp 

A dA 

The integration of /3 over dA is elementary and we arrive at 



(98) 



S(2,t)=p 3 (t). (99) 



This is the sought-after second-order differential equation. The coefficients are given 
in the equal mass case by l28l[l0l 

p (t) =t(t- to 2 ) (t - 9m 2 ) , p 2 (t) = t - 3m 2 , 

pi{t) = 3t 2 -20ttn 2 + 9tn 4 , p 3 (t) = -6n 2 . (100) 
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The coefficients for the unequal mass case can be found in 1321 . 



3.2 Differential equations 

The ideas of the previous subsection can be generalised to arbitrary Feynman in- 
tegrals. For a given Feynman integral let us pick one variable t from the set of the 
Lorentz invariant quantities (pj + p^) 2 and the internal masses squared m 2 . Let us 
write 



d»/V-(/+l)D/2 

<* = <»[ n*?" i V">/' ■ (101) 



Vf-l 

I 

\j =1 



The subscript t indicates that ft), depends on t through & . The Feynman integral is 
then simply 

I G = J (Or (102) 
A 

We seek an ordinary linear differential equation with respect to the variable t for the 
Feynman integral Iq. We start to look for a differential equation of the form 

L {r) (o,=dfi, (103) 

where 



is a Picard-Fuchs operator of order r. Suppose an equation of the form as in eq. ( 1103b 
exists. Following the same steps as in section ( 13. \\ we arrive at 

L {r] I G = J p. (105) 

dA 

The right-hand side corresponds to simpler Feynman integrals, where one propa- 
gator has been contracted. The coefficients of the Picard-Fuchs operator and the 
coefficients of the form /3 can be found by solving a linear system of equations | 
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